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Quasars as probes of the submillimetre cosmos at z > 5: 
I. Preliminary SCUBA photometry 
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ABSTRACT 

We present submillimetre (submm) continuum observations of a sample of some of 
the highest redshift quasars currently known, made with the SCUBA bolometer array 
on the JCMT. 12 of the sample have redshifts greater than five, and two have z > 
6; the median redshift of the sample is 5.3. Two of the z > 5 objects are strong 
(6a) detections, and are bright sources with Ss^o /j,m > lOmJy. Another firm (5a) 
detection is obtained for the z = 5.7 quasar SDSS J1044-0125; and SDSS J1306+0356, 
at z = 6.0, is detected with a signal-to-noise ratio «4. We have obtained sensitive 
(<7 «1.5mJy) upper limits for much of the remainder of the sample, including the 
z = 6.3 quasar SDSS J1030+0524. 

Submm spectral indices measured for two of the sources (a « 3.3) are consistent 
with thermal reradiation from dust, rather than from synchrotron emission. Sensi- 
tive upper limits at 450/im imply that the dust is cool, requiring large dust masses 
(1O 8_9 M0) to account for the observed fluxes — suggesting substantial prior star for- 
mation, even at z — 6 when the Universe was only l.OGyr old. 
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1 INTRODUCTION 

Forty years ago, an accurate optical identification was ob- 
tained for the curious radio source 3C273 (Hazard, Mackey 
& Shimmins, 1963). Subsequent spectroscopy revealed an 
unexpectedly high redshift, z — 0.16, heralding a dramatic 
new type of astrophysical phenomenon — dubbed the "quasi- 
stellar object" (QSO). Ever since, the extremely luminous, 
compact emission from quasars has rendered them almost 
unrivalled as beacons to the distant Universe — whether as 
spotlights shining through absorption systems, pinpoints fix- 
ing their host galaxies or as tracers of large scale structure. 

Recently, however, dedicated searches for galaxies at 
the highest redshifts have unveiled new populations of star- 
forming sources at z > 5 and z > 6 (e.g. Hu et al. 2002; 
Stanway, Bunker & McMahon 2003). Nevertheless, the lat- 
est AGN (Active Galactic Nuclei) surveys keep pace, and 
efficiently yield quasars at z > 5 (refs. under Fan et al; 
Sharp et al. 2001), with the promise of many more to be 
discovered over the coming years. That AGN remain impor- 
tant cosmological tools was recently highlighted by the first 



evidence for Gunn-Peterson absorption, discovered in the 
spectrum of a z = 6.3 quasar (Becker et al. 2001). The high- 
redshift quasars of this new generation stand tantalizingly 
as signposts to the epoch of reionization. 

The very existence of luminous quasars at z > 5, in- 
deed, poses a challenge to structure formation theory. Their 
supermassive black hole engines (M ~ 1O 9 M0), and the 
fuel reservoir required to supply them (M ~ 100 Mq yr _1 ), 
constitute a large mass concentration already in place only 
a gigayear after recombination. The Cold Dark Matter cos- 
mogony, in contrast, predicts that the first objects to form 
were of predominantly low mass (~ 1O 6 M ). High-redshift 
quasars therefore offer a a test of hierarchical cosmology 
operating at its extreme (Efstathiou & Rees, 1988; Turner, 
1991). The signatures of high overdensity, for example the 
biased evolution of the surrounding region, can in principle 
be detected observationally by targetting quasars and their 
environs (e.g. in the submm, Ivison et al. 2000). 

Here, we report continuum submillimetre observations, 
made with the SCUBA bolometer array on the James Clerk 
Maxwell Telescope, of a sample of some of the highest- 
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redshift quasars currently known. These observations form 
part of a larger project to establish the fiducial proper- 
ties of cool dust emission from radio-quiet, optically-selected 
quasars over a wide range in redshift (z < 1: Isaak et al., in 
prep.; z = 2: Priddey et al. 2003a; z > 4: McMahon et al. 
1999, Isaak et al. 2002). Our earlier work suggested that a 
significant fraction of optically-bright quasars are ultra- or 
even hyper-luminous at far-infrared wavelengths. Addition- 
ally, no significant difference between the z — 2 and z > 4 
samples was detected, contrasting with the SCUBA radio 
galaxy survey by Archibald et al. (2001). Several lines of 
evidence suggest that the submm light is powered by star 
formation, and multiband follow-up is in progress to test this 
hypothesis. Observation beyond z > 5, reaching yet further 
back in cosmic time, is a vital extension to the project. Arc 
conditions at the highest redshifts favourable, or adverse, 
to the formation of luminous, dusty sources? Do we catch 
the quasars during an even more turbulent phase of youth 
(gas-rich, violently star- forming), or does there exist a red- 
shift cut-off beyond which too little dust has formed for the 
object to be a SCUBA source? 

This Letter is the first in a series of papers discussing 
submm sources at z > 5. Here we present initial submil- 
limetre photometry of the sample, leaving a presentation of 
follow-up data, and a more thorough account of interpreta- 
tion, to forthcoming works (Isaak et al. in prep.; Robson et 
al. in prep.). Throughout, we assume a flat, A-dominated 
cosmology 57m=0.3, Q.a=0.7 with _ffo=65kms~ 1 Mpc -1 . 
Additional quantities (in parentheses) are for an Einstein-de 
Sitter (EdS) universe with .ffo=50kms _1 Mpc -1 . 



2 OBSERVATIONS 

The observed quasars are listed in Table 1. 

2.1 The sample 

The parent sample consisted of all quasars above a redshift 
of 4.90 that were known at the time of observation (i.e. late 
2001). The targets were drawn from the Sloan Digital Sky 
Survey (SDSS: see Fan et al. refs.), the Isaac Newton Wide 
Field Survey (WFS: Sharp et al., 2001) and the single object 
(RD J0301+0020) reported by Stern et al. (2000). Targets 
were prioritized solely on the basis of redshift, with target 
visibility providing the only additional constraint. 

Absolute B magnitudes are calculated from 1450A con- 
tinuum fluxes obtained from published spectra, assuming 
an optical spectral index a = —0.5. In Table 1, the primary 
values are for the A cosmology, and those for Q.m = 1 are 
given in parentheses. In our previous submm quasar surveys 
(McMahon et al., 1999; Isaak et al., 2002; Priddey et al., 
2003), we had targetted quasars brighter than Mb = —27.5 
(based on the EdS cosmology). At z > 5, too few sources 
are known for this stringent criterion to be maintained: the 
median and mean magnitude are each now Mb = —26.6 
(again EdS). 

2.2 JCMT-SCUBA submillimetre data 

Data for SDSS J1044— 0125 were obtained through service 
observing time in July 2000 and January 2001. The bulk 



of the rest of the data were obtained during the period au- 
tumn 2001-summer 2002, through a combination of "flexi- 
ble" scheduling, the scheduled time itself and Director's Dis- 
cretionary Time. 

SCUBA was employed in photometry mode, with the 
wide 850:450 filter set, and a standard 60 arcsec chop in 
azimuth at 7.8Hz. The source was placed on the central 
bolometer (H7, C14), and the median of the remaining 
(quiet) bolometers was used for additional sky removal. 
Flux calibration was performed against the planets Mars 
and Uranus, and against the standard continuum calibrators 
CRL618, OH231.8, IRC10216 and 16293-2422. Telescope 
pointing was checked frequently. Sky opacity was monitored 
via regular skydips at 850 and 450pm, and continuously at 
225GHz via the JCMT Water Vapour Monitor and the CSO 
Tau Meter: over the whole set of observations, T225 ranged 
between 0.04 < T225 < 0.14. Data were reduced with both 
the ORAC-DR pipeline and the SURF software. 

Both SDSS J1044-0125 and SDSS J1306+0356 were 
reobserved, after showing some inconsistencies in their ini- 
tial datasets: the detections were confirmed in each case. For 
example, during initial observation of SDSS J1306+0356, in 
November 2001, the final two blocks of integrations were in- 
consistent with the data taken over the preceding few hours. 
Their mean signals were zero or negative, where the others 
had throughout been strongly positive. This marked change 
in consistency coincided not only with deteriorating weather 
(seeing poorer than 4 arcsec), but with sunrise, and conse- 
quent fleeting changes in atmosphere, dish figure and ab- 
solute pointing accuracy. We reobserved the quasar under 
good, stable conditions (seeing 0.5 arcsec, t225gh z = 0.065), 
conditions in March 2002. The signal thereby obtained was 
consistent with the first seven blocks from the previous 
November. The final flux quoted in Table 1 therefore con- 
sists of the November data minus the last two blocks, plus 
the March data. Adding back the suspect data would give 
S , 850/jm=3.1±0.9mJy — a signal-to-noise which is still greater 
than three. 

SCUBA consists of two detector arrays, dedicated re- 
spectively to long- and short-wavelength observations — in 
this case 850 and 450/im. The atmospheric transmission de- 
creases with frequency, hence the long-wave array is the 
usual primary instrument. Nevertheless, 450/im data are ob- 
tained gratis, and they are here of particular interest, for 
they provide important constraints on the physical temper- 
ature of the dust. 450/im fluxes are thus included in Table 1 
for reference, however none of the sources was detected in the 
band. Observations were made under a wide range of atmo- 
spheric conditions (450/im opacities 0.75 < T4 5 o Mm < 3.0), 
resulting in a very heterogeneous set of RMS values, between 
10 and lOOmJy. We emphasise that accurate calibration of 
these short-wavelength observations is possible only under 
the very best observing conditions. For comparison, the rel- 
ative errors in flux calibration with the 850 and 450 arrays 
are typically 5-10 percent and ~ 20 percent, respectively. 



3 RESULTS 

Two of the targets, SDSS J0756+4104 and SDSS 
J0338+0021 (both 2 = 5.1) are strikingly bright submm 
sources, with Ssso/jm > lOmJy. The latter was also detected 
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Table 1. Quasars at z > 5 observed with SCUBA 
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Refs: 1. Fan et al. 2001; 2. Fan et al. 2000b; 3. Stern et al. 2000; 4. Anderson et al. 2001; 5. Zheng et al. 2000; 
6. Sharp et al. 2001; 7. Fan et al. 2000a; 8. Fan et al. 1999 

a Observations made during high atmospheric opacity (Weather Grade 3): short-wave data consequently very poor 
6 Radio source: Si.4GHz=17.7mJy (NVSS catalogue) 



at 1.2mm by Carilli et al. (2001), Si. 2mm = 3.7mJy. The 
850pim flux is consistent with a thermal spectrum (see Fig- 
ure 1 and Section 4.1). SDSS J1044-0125 (z = 5.7) is a 
moderately bright detection, Ssso = 7mJy. Its 1.2mm flux 
of around 2mJy (Isaak et al., in prep.) is consistent with 
the steep Rayleigh-Jeans tail of a thermal greybody spec- 
trum from dust. This source and its environs have formed 
the target of a range of follow-up observation, reported in a 
companion paper (Isaak et al., in prep.). It is notable that 
SDSS1044 is a broad absorption line (BAL) quasar, a fact 
drawn upon to account for its X-ray weakness, relative to 
the optical, as measured by Brandt et al. (2001). It is con- 
ceivable that the gaseous outflow responsible for the optical 
absorption also absorbs the X-rays; and that dust embed- 
ded within the gas gives rise to the submm emission. SDSS 
J 1306+0536 (z = 6.0) was detected at 850^m with a sig- 
nificance of 3.7a (see also Section 2.2). 

RD J0301+0020 (z=5.50) was detected at 1.2mm in a 
very deep IRAM-MAMBO observation by Bertoldi & Cox 
(2002): the 850/im limit presented here is consistent with 
a thermal spectrum given their Si. 2 = 0.87 ± 0.20mJy. RD 
J0301+0020 stands out from the other z > 5 quasars by 
virtue of its extremely low optical luminosity, Mb = —23.1. 
Using the Elvis et al. (1994) bolometric correction from vLb 
(Cb = 12) implies Lboi ~ 2 x 10 12 L . This is comparable to 
the far-infrared luminosity derived from the millimetric flux, 
if we assume a canonical cool dust spectrum (T = 40 — 50K). 
Assuming that the quasar (and not some foreground object) 
is the MAMBO source, then either the AGN is absorbed, and 
we are seeing reprocessed emission in the far-infrared, or the 
Lfir derives from star formation (or, indeed, a combination 
of the two). 

3.1 Caveat emptor: the effects of lensing 

In neither Table 1 nor Table 2 do we attempt to correct the 
observed fluxes, and their derived quantities, for a gravita- 



tional lensing magnification. At z > 5, one would expect the 
optical depth to lensing to be significant; Wyithe & Loeb 
(2002) additionally point out that the z > 5 SDSS survey 
selects quasars on the steep, bright tail of the luminosity 
function, rendering it prone to magnification bias. Gravita- 
tional lensing is therefore a serious problem whose effects 
must be addressed on a source-by-source basis. For now we 
can only do so statistically: for our two z > 5.5 detections, 
SDSS J1306 and SDSS J1044, Wyithe & Loeb (2002) esti- 
mate a 7-30 percent probability of multiple imaging, with 
median (mean) magnifications between 1.1 (5) and 1.2 (30). 

Schwartz (2002) examined the Chandra- ACIS images 
of the three highest-redshift objects in this sample. The six 
photons comprising the SDSS J1030 detection are widely 
distributed, statistically inconsistent with a point source. 
This is a possible signature of lensing, however the counts are 
too few to differentiate between a single extended source or 
multiple point sources. In contrast, Fan et al. (2003) report 
that HST-ACS images of SDSS J1030 and J0836 are con- 
sistent with unresolved point sources. Their ground-based 
imaging of SDSS J1306 promotes a similar conclusion. 



4 DISCUSSION 

4.1 Dust, metals and star formation at z > 5 

The mm-submm spectral indices measured for SDSS 
J1044-0125 and J0338+0021— each a™^ mm « 3.3— are 
consistent with the Rayleigh-Jeans tail of thermal emis- 
sion from cool dust. Such a spectral energy distribution 
(SED) would be expected to peak at a rest wavelength 
« 50— 100/im; similar, for example, to star- forming galaxies 
such as M82. Indeed, existing 450 and 350pm detections of 
z ~ 4 quasars suggest that the SED tends to a plateau (see 
e.g. Priddey & McMahon, 2001 (PM01)) at around lOO^m 
rest-frame: all but the brightest sources would lie below 
the short-wavelength limit that we could detect. At z > 5, 
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Table 2. Derived properties of the sample (detections and deep (a < 1.5mjy) limits). Numbers in parentheses are for the EdS cosmology 
(except for M*(min), where the numbers for each cosmology are approximately equal). RD0301 is counted among the deep 850/im 
non-detections, but we have used its 1.2mm flux to calculate the parameters. N.B. Lfir and M* are not independent. 

Source z t(oo) - t(z) M d M*(min) L FIH M bh M acc 

Gyr 1O 8 M Mgyr 1 10 10 L Q 10 9 M Meyr 1 
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Figure 1. Constraining the dust temperature of SDSS0338. The 
1.2mm point is from Carilli et al. (2001), the 850 and 450/im 
from this work. Curves are for the maximum dust temperatures 
(labelled) consistent with the 450/im upper limit (plotted as signal 
plus 2-a with a c-length bar), for three greybody indices: (3 = 2.0 
(solid), (3 = 1.5 (dashed) and /3 = 1.0 (dotted). The steep submm- 
mm slope renders a better fit for a high f3, thus favouring a low 
temperature. 



then, a 450/itm detection would imply that the dust is hot; 
conversely, the deep limits obtained for the bright 850/im 
detections are valuable constraints on their dust temper- 
ature. Again considering J0338+0021, 450 upper limit re- 
quires T < 100K for a greybody index 13=1, and T < 50K 
for f3 = 2, the higher j3 favoured by the steep mm-submm 
slope (Figure 1) . (However we warn of the potentially consid- 
erable systematic uncertainties in the 450/im calibration.) In 
the following, we shall therefore assume the SED of PM01, 
namely (3=2, T=A0K. 

Temperature uncertainty notwithstanding, determina- 
tion of the dust opacity (k) has, in the past, presented ad- 
ditional uncertainties to a derivation of a dust mass from 
a submm flux. We adopt a value ft(125/im)=30 cm 2 g _1 , at 
a normalisation wavelength (125pm) corresponding to ob- 
served 850/im at z=5.8. This is consistent with the k(125) 
determined by Hildebrand (1983). Extrapolated to longer 



wavelengths assuming k oc A~ 2 , it is also consistent with 
a Galactic mixture of Draine & Lee (1984) grains as well 
as with more recent measurements, e.g. the 850/im value of 
James et al. (2002). For our z > 5 sample, it yields masses 
in the range 10 8 " 9 M Q (Table 2). 

The youth of the universe — lGyr at z = 5 — imposes 
a constraint on the mechanisms of dust production, for it 
is comparable to the evolutionary timescale of the stars in 
whose atmospheres the dust is believed to condense. An al- 
ternate method of manufacture, on the winds of Type II 
supernovae, requires much shorter timescales, but its effi- 
cacy has yet to be proven. Crudely assuming an overall dust 
enrichment rate per mass forming stars of around 1 per- 
cent, and an absolute upper limit of lGyr available for star 
formation, then a dust mass of 10 8-9 M would require a 
minimum average star formation rate of ~ lOOM yr _1 . 
It is thus plausible that star-forming activity is responsible 
for some or all of the observed submm luminosity. In Ta- 
ble 2, the quantity M* (min) is the minimum star formation 
rate necessary to produce the estimated dust masses, aver- 
aged the lifetime of the universe. Alternatively, Lpir could 
be used to estimate the instantaneous star formation rate, 

; ; M * — r = * ~R , where 'I' is a constant of order unity, 

M yr _1 io 10 L© J ' 

depending on the stellar initial mass function (IMF). This 
quantity represents an upper limit on the star formation 
rate, due to the potential of AGN contribution to dust heat- 
ing. Nevertheless it is plausible that the dust we observe was 
formed during production of a substantial fraction of stars 
in the quasar's host galaxy. 

Morgan & Edmunds (2003) present a detailed model 
of dust synthesis and discuss its implications for very high 
redshift submm sources such as the current sample. Their 
model implies that, in the absence of supernovae as a source 
of dust, the star formation efficiency must be very high to 
yield, by z = 5, the masses we observe. 



4.2 The growth of supermassive black holes 

The four highest-redshift quasars in the sample are all very 
luminous in the optical, Mb ~ —28. The assumption that 
this luminosity is supplied by Eddington-limited accretion 
then requires a black hole of at least 5xlO 9 M (Table 2). 
How does a black hole of such mass form? Imposing the max- 
imum available time of ~lGyr implies an average accretion 
rate of 5 M yr _1 . In comparison, assuming an (optimistic) 
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radiative efficiency e = 0.1 (where Lboi = tMc 2 ), an (instan- 
taneous) accretion rate of 100 Mq yr -1 is required to fuel 
the bolometric power inferred from the optical. The e-folding 
time for Eddington accretion is w 0.5eGyr; starting from a 
10 3 ~ 6 M seed (e.g. Haehnelt, Natarajan & Rees, 1998), an 
observed redshift of z~6 would require a formation redshift 
2=10-18 in the A cosmology, but an upper limit longer than 
the age of the Universe in the matter-dominated cosmology. 

In the local universe, quiescent supermassive black holes 
and their surrounding stellar bulges correlate in their mass 
and velocity dispersion (Magorrian et al. 1998; Ferrarese 
& Merrett 2000; Gebhardt et al. 2000). It is likely, there- 
fore, that their formation processes — the one through accre- 
tion and an AGN phase, the other perhaps through a dusty 
starburst — are closely linked. What can we conclude con- 
cerning the coevolution of a quasar and its host galaxy, given 
these data? Both PM01 and Archibald et al. (2002) consider 
this matter. They propose a simple model in which a black 
hole exponentiating according to the Eddington rate, and 
an elliptical host galaxy forming stars and dust at roughly 
constant rate, conspire to produce a submm-luminous AGN 
just before the gaseous fuel supply is exhausted. In a future 
paper, we will consider a wider range of possibilities. 



5 CONCLUSIONS AND FUTURE WORK 

We have presented SCUBA 850/im observations a sample 
of quasars at z > 5 — where the youth of the Universe it- 
self («lGyr) starts to become a constraint on the histo- 
ries of accretion and star formation — detecting four sources 
at 850/im. The few 1.2mm-850/im spectral indices that we 
can determine are consistent with thermal reradiation from 
dust; and deep 450/um limits imply that the dust is cool 
(T < 100K), consistent with their siblings at z = 4 (Priddey 

6 McMahon, 2001). Although the sample is much smaller, 
and its median optical luminosity fainter, its submm proper- 
ties seem tentatively similar to those of radio-quiet quasars 
we have studied at z — 4 and z = 2 (Isaak et al. 2002; 
Priddey et al. 2003a). 

It is plausible that we are observing extreme, high-u 
peaks in the overdensity distribution. The quasar host galax- 
ies are likely massive, gas-rich and actively forming stars. 
We are pursuing a range of follow-up observations designed 
to investigate such a scenario. For example, submillimetre 
imaging of the fields of high-redshift AGN has been used 
(e.g. Ivison et al. 2000) to test for biased galaxy formation. 
In this spirit, we have obtained SCUBA jiggle maps of three 
of the sources from the current sample, and our findings will 
be reported in a forthcoming paper (Isaak et al. in prep.). 

An important quantifier of the evolutionary state of 
the sources themselves will be obtained through detection 
of carbon monoxide emission lines, which would reveal the 
existence of reservoirs of molecular gas indicative of star for- 
mation. Simultaneously, it will be essential to improve the 
constraints on the infrared/submm SEDs of these quasars, 
to confirm the presence of dust, differentiate between AGN- 
and starburst-powered components and improve estimates 
of the dust mass and star formation rate. And of no less im- 
portance is the need to obtain high resolution optical/near- 
infrared images of the quasars, to assess the probable effects 
of gravitational lensing. 



It is an oft-cited consequence of the negative submm K- 
correction that it hypothetically permits an unidentified 
SCUBA source to lie at a redshift as high as ten, without 
being intrinsically more luminous than a local object like 
Arp220. The difficulty of optically identifying such a source 
has prevented this claim becoming more than hypothetical. 
Now, however, the submillimetre study of samples of quasars 
at z > 5 provides a ready means of inferring the properties- 
indeed the very existence- of dust within one gigayear of the 
Big Bang. 
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